Stellar Structure and Evolution Stellar Models 


Due February 11 


Your task is to compute a stellar model for a chemically homogeneous star using the Matching 
Point Method. You may choose whatever mass you wish, though it may be easier to get 
your model to converge if M > 2 Mo. Assume that the composition of the star is X = 0.75, 
Y = 0.23, Z = 0.02, and use at least Z 50 shells. You may make the following assumptions: 


a) Neglect any degeneracy. Except for the very lowest mass stars, degeneracy is not important 
on the main sequence. The equation of state will therefore contain just two terms: one for 
an ideal gas (both ions and electrons) and one for radiation pressure. 


b) Neglect all energy sources except for the proton-proton chain and the CNO bi-cycle. The 


energy generated by these reactions (in ergs s7! cm™?) can be approximated via 
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where Ty is the gas temperature in units of billions of degrees. 
c) For opacity, you may interpolate in the tables given by 


http://cdsweb.u-strasbg.fr/topbase/The0P. html 
or 
http://opalopacity.1llnl.gov 


or use the rough approximations given in the notes. If you choose the latter, you can neglect 
H` opacity (for stars more massive than about 1.5 Mg), and use 
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where the opacity is given in terms of cm? g~!. 

Your output should consist (at minimum) of plots of pressure, temperature, luminosity and 
radius versus the interior mass of the star. Optionally, you may also make plots of K, €pp, 
ECNO; V, Vrad: Vad; and/or compare any variable against any other variable. 


